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I. INTRODUCTION 

Cosmological inflation [1] has become the dominant paradigm to understand the initial conditions for Cosmic 
Microwave Background (CMB) anisotropies and structure formation. In the inflationary picture, the primordial 
cosmological perturbations are created from quantum fluctuations “redshifted” out of the horizon during an early 
period of accelerated expansion of the universe, where they remain “frozen”. They are observable as temperature 
anisotropies in the CMB. They were first detected by the COsmic Background Explorer (COBE) satellite [2, 3, 
4]and then with spectacular accuracy most recently by the the Wilkinson Microwave Anisotropy Probe (WMAP) 
mission [5, 6]. Since the observed cosmological perturbations are of the order of 10“^, one might think that first-order 
perturbation theory will be adequate for all comparison with observations. That may not be the case however, as 
the Planck satellite [7] and its successors may be sensitive to the non-Gaussianity of the cosmological perturbations 
at the level of second- or higher-order perturbation theory [8]. Statistics like the bispectrum and the trispectrum of 
the CMB can then be used to assess the level of primordial non-Gaussianity on various cosmological scales and to 
discriminate it from the one induced by secondary anisotropies and systematic effects [8, 9, 10, 11]. Therefore, it is 
of fundamental importance to provide accurate theoretical predictions for the statistics of CMB anisotropies as left 
imprinted by the primordial seeds originated during or immediately after inflation. Steps towards this goal have been 
taken in Refs. [12, 13, 14, 15, 16, 17] at the level of second order perturbation theory. 

In this paper we provide the first computation of the complete second-order radiation transfer function for CMB 
anisotropies on large angular scales in a flat A Cold Dark Matter (ACDM) model. Our calculation includes the 
second-order generalization of the Sachs-Wolfe as well as of the (both early and late) Integrated Sachs-Wolfe effects 
and is valid for a generic set of initial conditions specifying the level of primordial non-Gaussianity. It also accounts 
for the contribution of second-order tensor perturbations. Our results which are complete as fas as large angular scales 
are concerned, represent the first step towards the computation of the full second-order radiation transfer function 
for GMB anisotropies. Its computation also on small angular scales, which involves several non-linear effects, such 
as gravitational leasing, Shapiro time-delay, Rees-Sciama effects, etc. and, of course, dealing with the non-linear 
Boltzmann equation for the photon-matter fluid (some specific aspects of which have been studied in Ref. [18]), will 
be the subject of a future publication [19]. 

The general expression for CMB anisotropies up to second order from gravitational perturbations has been obtained 
in Refs. [12, 13]. At large angular scales the extension up to second order of the Sachs-Wolfe effect [20] has been 
derived in the Poisson gauge (for adiabatic perturbations) in Refs. [8, 21] (and for a gauge-invariant generalization at 
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second-order see Ref. [14]) ^ 
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(1.1) 


where (/)* = - I - ^ 2 */^ =(/?*-!- 02*/2 is the gravitational potential at emission on the last scattering surface 

and is split into linear and second-order parts (see Eq. (2.1)). Such an expression allows to single out the initial 
(at last scattering) non-linearities of the large-scale CMB anisotropies, including the primordial non-Gaussianity of 
cosmological perturbations, as it has been computed in Refs. [8, 14, 21, 22]. In Eq. (1.1) the large-scale second-order 
gravitational potential (j) 2 * at emission contains the non-linearity arising from inflation and from the post-inflationary 
epoch yielding [8, 14, 21, 22] 

at 1 1 2 ^ 2 

where 

, ( 1 - 3 ) 


1C = iov-4a,aj(5v*a,<^*) - v-^ ( H 


and Oni is a parameter specifying the level of primordial non-Gaussianity which depends upon the particular scenario 
for the generation of cosmological perturbations, as we will discuss later in more detail. For example, for standard 
single-field slow-roll models of inflation Oni = 1-1- 0(e,r]) [23, 24, 25], where e and 77 are the standard slow-roll 
parameters (see e.g. Ref. [1]). 

In this work we will focus on another contribution to the large-scale CMB anisotropies which is due to the evolution 
of the gravitational potentials from the epoch of last scattering up to now and it is encoded in a number of integrated 
terms [13] 

AT d 

— (n)=y dT]—{(l) + ^p){x,r]) 

where (j) = (jji + ^ 2/2 and 'tp = ip i + '02/2. In Eq. (1.4) n denotes the direction of photon momentum and rj is 
the conformal time, rjQ and ry* being the present and the last scattering times, respectively. Note that the partial 
time derivative of the perturbations must be evaluated along the line of sight. These integrated terms represent the 
gravitational redshift due to the time variation of the metric during the photon travel from the last scattering surface 
to the observer. The linear contribution is the so-called Integrated Sachs-Wolfe (ISW) effect and in a flat matter 
dominated (Einstein-de Sitter) Universe it vanishes due to the constancy in time of the gravitational potentials. 
However, in the presence of a cosmological constant (or more general dark energy component) it actually gives a 
large-scale contribution due to the decay in time of the potentials: this is the so-called late ISW effect. The early ISW 
effect similarly arises because of the residual dynamical role of the radiation component at times close to the last¬ 
scattering epoch; this also makes the gravitational potentials evolve in time, leading to a non-negligible contribution 
to the large-scale anisotropy (see e.g. Ref. [26]). The second-order generalization of the ISW effect is characterized 
by some pecnliar properties. Even in a flat matter-dominated Universe, when the linear counterpart vanishes because 
the gravitational potential is constant, the evolution of the perturbations in the non-linear regime generates a non¬ 
vanishing contribution to the GMB anisotropies, the so-called Rees-Sciama effect [27]. In this paper we will compute 
the second-order ISW generated in a universe with a non-vanishing cosmological constant term. Unlike the usual 
Rees-Sciama effect, a contribution to the CMB anisotropies on large scales is generated, and therefore for the rest of 
the paper by second-order ISW we will refer to the large scale part of Eq. (1.4). Such an effect has not been considered 
so far, and thus a major part of this paper will deal with its study. An important feature is that it is related to possible 
primordial non-Gaussianity in the energy-density perturbations produced in the early universe. In a separate paper 
we will give a quantitative study of its signatures in the bispectrum of the CMB anisotropies [19]. Notice that in 
Eq. (1.4) we have also indicated additional contributions to CMB anisotropies which arise from integrated terms 
quadratic in the linear metric perturbations [12, 13]. We will discuss and take them into account later. Finally in 
writing Eq. (1.4) we have dropped possible gravitational redshifts from second-order vector modes which are expected 
to give a negligible contribution on large scales [13]. 

This paper serves for different purposes. First it contains an independent part about the evolution of the second- 
oder perturbations in a flat universe filled with matter and a cosmological constant term. Second it discusses the 


x= —rifrin—rjl 


(first — order) ^ 


(1.4) 


^ For a fully non-linear generalization of the Sachs-Wolfe effect see Ref. [22]. 
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generation of CMB anisotropies due to the second-order ISW and its connection to primordial non-Gaussianity. The 
final achievement is that of providing the full second-order radiation transfer function for CMB anisotropies on large 
angular scales from the second-order Sachs-Wolfe and ISW effects. Our formalism is the starting point of an accurate 
evaluation of different effects on the CMB anisotropies, including the exact expression for the CMB bispectrum on 
large scales [19]. 

The paper is organized as follows. In Section II we present the main equations and solutions for the evolution of 
the gravitational potentials in a ACDM universe, accounting for primordial non-Gaussianity in the initial conditions 
(see also Appendix A). In Section III we derive the expression for the second-order ISW by taking the large-scale 
contributions to Eq. (1.4) and we derive an expression of the CMB anisotropies expanded in multipoles (with details 
in Appendix B). In Section IV we study the contribution to the CMB anisotropies from second-order tensor modes. 
Finally in Section V we provide the full expression of the second-order radiation transfer function for large-scale CMB 
anisotropies. In this section we summarize our final formulae as a list and we provide the reader with all the necessary 
definitions. Section V contains also our conclusions. 

II. SECOND-ORDER PERTURBATIONS OF A FLAT ACDM UNIVERSE IN THE POISSON GAUGE 

A. Evolution of the second-order gravitational potentials 

We consider a spatially flat Universe filled with a cosmological constant A and a non-relativistic pressureless fluid 
of Cold Dark Matter (CDM) plus baryons, whose energy-momentum tensor reads = pu^u^. The radiation 
component, which also plays a role during the early stages following the last-scattering epoch, will be considered in 
Section IIIC, in connection with the early ISW effect. We will work in the so-called Poisson gauge which, to linear 
order, reduces to the Newtonian gauge. Using the conformal time dp = dt/a, the perturbed line element around a 
spatially flat FRW background reads 

ds^ = a^{p){—{l + 2(j))dp^ + 2ojidpdx'‘ + [(1 — 2'>p)6ij + Xij]dx^dx^} . (2.1) 

Here each perturbation quantity is expanded into a first-order (linear) part and a second-order contribution in the 
same way as for the gravitational potential (j) = (f>i + 02/2. In the Poisson gauge the shift perturbation uji is a pure 
vector, i.e. diUt^' = 0, while Xij is a tensor mode {i.e. divergence-free and traceless d^'Xij = 0, x^'i = 0). la this 
paper we are mainly interested in the evolution of the second-order scalar potentials. Moreover we will neglect linear 
vector modes since they are not produced in standard mechanisms for the generation of cosmological perturbations (as 
inflation), and we also neglect tensor modes at linear order, since they give a negligible contribution to the bispectrum 
arising from the integrated effects (1.4). In Section IV and Appendix A, we will report also the computation for the 
second-order vector and tensor modes following the results of Ref. [28]. As for the matter component we split the 
mass density into a homogeneous p{p) and a perturbed part as p{x,p) = p{p){l -f di -I- 02 / 2 ) and we write the four 
velocity as = (Sq + v^)/a with = —1 and + V 2 / 2 . 

The Friedmann background equations are 370^ = (87rGp -I- A) and p' = —ZTip, where a prime stands for differen¬ 
tiation with respect to conformal time, and Ti. = a'/a. 

Let us briefly recall the results for the linear perturbations in the case of a non-vanishing cosmological A term. At 
linear order the traceless part of the (i-j)-components of Einstein equations gives (pi = tpi = tp. Its trace gives the 
evolution equation for the linear scalar potential p 

p" + mp' -b a^kp = 0 . (2.2) 

Selecting only the growing mode solution one can write 

V 3 (x,r 7 ) = 5 ( 77 ) 1 ^ 0 (x), (2.3) 

where pn is the peculiar gravitational potential linearly extrapolated to the present time (po) and g{p) = D+{p)/a{p) 
is the so called growth-suppression factor, where D^{p) is the linear growing-mode of density fluctuations in the 
Newtonian limit. The exact form of g can be found in Refs. [29, 30, 31]. In the A = 0 case g = 1. A very good 
approximation for g as a function of redshift z is given in Refs. [29, 30] 

5 cx - Ha + (1 + H„/2) (1 + Ha/70)] , (2.4) 

with Hm = Hom(l + z)^/E‘^{z), Ha = ^oa/E‘^{z), E{z) = (1 + z)H{z)jHo = [Hom(l + ^)^ + Hqa] and Hq™, 
Hqa = 1 ~ Hom, the present-day density parameters of non-relativistic matter and cosmological constant, respectively. 
We will normalize the growth-suppression factor so that g{z = 0) = 1. 
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Let us now consider the second-order scalar perturbations. For simplicity, we adopt the notation 


02 = $, (2.5) 

^^2 = «'. (2.6) 


The evolution equation for the second-order gravitational potential 4^ is obtained from the trace of the (i-j)-Einstein 
equations: its detailed derivation will be given in Appendix A; here we just report it 


4'" -K 37f4'' -f a^A4' = S{r]), 


where S{ri) is the source term 


S{r]) = g^VLmH^ 




ifQ + 2 2 




12 ri 


+ 3 X V ^ (idVo'^i'Po) - 3V [d'^ipodjipo) 


^ ^ - (aVo^iV^o) 


and we have introduced the function 


f! ^ _ dlni:»+ , g'{r]) 


(2.7) 

( 2 . 8 ) 

(2.9) 


which can be written as a function of flm as /(flm) ~ [29, 30]. In Eq. (2.8) stands for the inverse of 

the Laplacian operator. 

The solution of Eq. (2.7) is then given by 


4'(r7) = 

9m 






( 2 . 10 ) 


Here 'i'+{ri) = g{r]) and ^'-(r/) = H{ri)/a‘^{ri) are the solutions of the homogeneous equation, W is the Wronskian 
which explicitly reads IF(? 7 ) = Ho/a^ with Wq = Ti-Qi/o + (3/2)Hom)) where the suffix ‘0’ stands for the value of the 
corresponding quantities at the present time. In particular notice that = '^{'gm) represents the initial condition 
taken deep in the matter dominated era on super-horizon scales, rjm being the epoch when full matter domination 
starts. It is such an initial value that must be properly determined in order to account for the primordial non- 
Gaussianity in the cosmological perturbations. The evolution of the gravitational potential $ is then obtained from 
the relation between 4' and 4) 


V2v24-= 


V2(5*(po5Vo) - 


( 2 . 11 ) 


which follows from the traceless part of the (t-j)-component of Einstein equations (see Appendix A). In particular 
deep in the matter dominated epoch Eq. (2.11) reduces to 

- Agly^V\l - ^gl, [v2(5Vo5*<^o) - , (2.12) 

where we have used that for ^ 1 f{g) 1, and where ipm = gm^Po is the value of the gravitational potential 
during matter domination, when the cosmological constant was still negligible, with g^ = g{Vm)-'^ Note that in this 
way Eq. (2.12) reduces to the relation for the two gravitational potentials 4> and d' first obtained in Ref. [25] for a 
matter-dominated Universe. 


B. Initial conditions from primordial non-Gaussianity 

Let us now discuss the key issue of the initial conditions which we conveniently fix at the time when the cosmological 
perturbations relevant today for Large Scale Structures and CMB anisotropies are well outside the Hubble radius, i.e. 


2 A good approximation, g™ 
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when the (comoving) wavelength A of a given perturbation mode is such that A ^ = a! ja being the horizon 

size. 

In the standard single-field inflationary model, the first seeds of density fluctuations are generated on super-horizon 
scales from the fluctuations of a scalar field, the inflaton [1]. Recently many other scenarios have been proposed 
as alternative mechanisms to generate such primordial seeds. These include, for example, the curvaton [32] and 
the inhomogeneous reheating scenarios [33], where essentially the first density fluctuations are produced through 
the fluctuations of a scalar field other than the inflaton. In order to follow the evolution on super-horizon scales of 
the density fluctuations coming from the various mechanisms, we use the curvature perturbation of uniform density 
hypersurfaces C = Ci + C 2/2 -b • • •, where Ci = — V'l ~ 'H5pi/p' and the expression for ^2 is given by [34] 

C2 = ~V'2 ~ ^ + ^C2 ) (2-13) 

P 


with 


AC2 = m 


hp' 5xp 
P' P' 


+ 2—+ 27i'!/’i) 
P 



(2.14) 


The crucial point is that the gauge-invariant curvature perturbation C remains constant on super-horizon scales after it 
has been generated during a primordial epoch and possible isocurvature perturbations are no longer present. Therefore, 
we may set the initial conditions at the time when C becomes constant. In particular, C 2 provides the necessary 
information about the “primordial” level of non-Gaussianity generated either during inflation, as in the standard 
scenario, or immediately after it, as in the curvaton scenario. Different scenarios are characterized by different values 
of C 2 . For example, in the standard single-field inflationary model C 2 = 2 (Ci) + O {e,r]) [23, 25], where e and p are 
the standard slow-roll parameters [1]. In general, we may parametrize the primordial non-Gaussianity level in terms 
of the conserved curvature perturbation as in Ref. [14] 

C 2 = 2a„i(Ci)" , (2.15) 


where the parameter Oni depends on the physics of a given scenario. For example in the standard scenario Oni — 1, 
while in the curvaton case Oni = (3/4r) — r/2, where r « {P(tIp)d is the relative curvaton contribution to the 
total energy density at curvaton decay [8, 35]. In the minimal picture for the inhomogeneous reheating scenario, 
Oni = 1/4. For other scenarios we refer the reader to Ref. [8]. One of the best tools to detect or constrain the 
primordial large-scale non-Gaussianity is through the analysis of the CMB anisotropies, for example by studying the 
bispectrum [8]. In that case the standard procedure is to introduce the non-linearity parameter f^i characterizing 
non-Gaussianity in the large-scale temperature anisotropies [8, 36, 37]. To give the feeling of the resulting size of /ni 
when jonij ^ 1, /ni — 5ani/3 (see Refs. [8, 14]). The conserved value of the curvature perturbation ( allows to set 
the initial conditions for the metric and matter perturbations accounting for the primordial contributions. At linear 
order during the matter-dominated epoch and on large scales Ci = —5(/?m/3, where ipm = g-m^o is the value of the 
gravitational potential during matter domination, when the cosmological constant was still negligible. Thus we can 
write 


OU 2 2 2 /r\ -1 r>\ 

C 2 = ^rn ~ 9m^0 * (2.16) 

On the other hand by using the expression for C 2 during the matter-dominated epoch together with Eq. (2.12) and 
the second-order (O-O)-component of Einstein equations (evaluated for a matter-dominated epoch), one can express 
in terms of (2 as (see Refs. [8, 14, 21, 25]) 

3 16 

= -gC2 + yT’m + TsJ~'^{d\mdiPm) “ {d\mdjPm) ■ (2.17) 

Therefore using Eq. (2.16) we find, for generic non-Gaussian initial conditions parametrized by ^2 = 2ani (Ci)^; 


= 2g“ 


= 2 g^ 


-^(Oni - 1) + 7^0 + ^ ^ (9Vo9i(/?o) - 3V '^did^ {d'^podjPo) 


-q(®ni “ 1) “ 1 - q ^ ((9Vo'9i7’o) - 3V ((9Vo'9i7’o)) 


(2.18) 


(2.19) 
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where Oni will always signal the presence of primordial non-Gaussianity according to our parametrization. The last 
of these equations (2.19) has been obtained from Eq. (2.12). 

After having determined the initial conditions from Eq. (2.10) we finally get 


d'(77) = 2g{r))gr, 


(-^(oni - 1) - 1^ (/?o “ ^ ^ {^"(po^^(po) - 3V (SVo^j^o)) 


+ [ fo + 


^Om 


Vm 


a- 


Kv) 


dfia^{fi)g{g)S{fi) 


' 'HTn 


( 2 . 20 ) 


f 


$(r;) = ^-( 77 )+ V(r?)(^o + + o ) [^ ^ {^"ipo^^^po) - 3V (d^ipodjipo)] . (2.21) 


Notice that we can rewrite the solutions for dt and $ in a compact form as 
'^(v) = (Bi{g)-2g{g)g^-^{ani-l)g{v)gm)(pl+(B2{g)-^g{g)g.. 


$( 77 ) = 


+ Bz{g)V ^9i9^(5Vo37V3o) + B4(?7)9Vo3iV3o , 

Bliv) + V(?7) - ‘2g{v)gm - y («nl - l)ff(?7)5m^ v^o + 

V"^ (aVo3i(/?o) - {d"(fodj(po) 


V ^ {d^ipodi(fo) - 3V (^Vo'^i'^o) 

( 2 . 22 ) 


B2iv) + o9^id) e(? 7 ) + - - -9{v)9 


+ B3{r])V ^did\d^ipodjipo) + B4(77)aVo3i(/?o , 


(2.23) 


where we have introduced Bi{g) = Hg ^ (/o + 30om/2) ^ 


Bi(ji) with the following definitions 


Biiv) 

Bsiv) 


~\ £7 _o / j~u2i 


/ df,H\fi){fiv)-iyCig,fj), B2iv)=2 6177^^(77) 2(7(77) - 1)^ - 3 + 30^(7)) ^(77,77) ,(2.24) 

d?7 fe(7 ?) + I) ^(t?,?)), 54(77 )=-/ dfiC{g,fi), (2.25) 


'r]m 


and 




C{g,g) =9 (j7)a(j7)|5(?7)W(?7) -5(»?)y^W(»7) 


(2.26) 


with 6 ( 77 ) = f{g)/ilm{v)- 

Let us comment on the results found here, in order to proceed further by computing the radiation transfer 
function for the second-order CMB anisotropies on large scales. In the expression for the second-order gravita¬ 
tional potentials (2.22) and (2.23) we recognize two contributions. The term which dominates on small scales, 
[B 3 {g)\/~‘^did^{d^ipodj(fo) + B 4 {g)d^(podi(po], which gives rise to the second-order Newtonian piece and is insensitive 
to any non-linearity in the initial conditions. In fact, for a vanishing cosmological constant A ^ 0, it is just this 
term which is responsible for the Rees-Sciama effect due to the non-linear evolution of the gravitational potentials. 
For the Einstein-de Sitter case, the signatures of the Newtonian Rees-Sciama effect in terms of the CMB bispectrum 
have been widely discussed in the literature [38, 39]. An analysis of the CMB bispectrum induced by the Newtonian 
Rees-Sciama effect in the presence of a non-vanishing cosmological constant A 0 has been given in Refs. [40, 41]. 

The remaining pieces in Eqs. (2.22) and (2.23) represent the novelty of our work. They correspond to contributions 
which tend to dominate on large scales with respect to those characterizing the Newtonian contribution, and whose 
origin is purely relativistic. The contribution of these terms to the CMB anisotropies through Eq. (1.4) will be mainly 
on large scales through the time dependence of the growth-suppression factor g{g) and it corresponds to the late ISW 
effect at second order. Moreover these are the pieces carrying the information on primordial non-Gaussianity, since 
they properly take into account in the initial conditions the primordial non-Gaussianity generated during inflation, as 
in the standard scenario, or immediately after it, as in the curvaton scenario. Terms of this type have been obtained 
for a flat matter-dominated (Einstein-de Sitter) Universe in Refs. [13, 42], but in this case there is no large-scale 
contribution through the integrated effect (1.4), because the different terms are constant in time, since 7 /( 77 ) = 1 
and Bi{g) = B 2 {g) = 0. In Ref. [15] second-order cosmological perturbations have been computed in the A 0 
case from the synchronous to the Poisson gauge, thus extending the analysis of Ref. [42], and the CMB temperature 
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anisotropies induced by metric perturbations have been also considered by applying the expressions of Ref. [13]. 
However, an important point to notice is that both Refs. [13, 42[ and Ref. [15] disregard any primordial non-linear 
contribution from inflation. ^ 


III. SECOND-ORDER ISW EFFECT IN A ACDM UNIVERSE: PRIMORDIAL NON-GAUSSIANITY 
AND SECOND-ORDER RADIATION TRANSFER FUNCTION 

A. Temperature anisotropies 

In order to give an expression for the second-order Integrated Sachs-Wolfe effect we will take the large-scale contribu¬ 
tions to Eq. (1.4). Thus in Eqs. (2.22),(2.23) for the second-order gravitational potentials we will drop the Newtonian 
terms that have already been examined in Refs. [40, 41] in connection to the Rees-Sciama effect. Accordingly we will 
take the second-order corrections (first — order)^ - quadratic in the linear perturbations - appearing in Eq. (1.4) from 
Refs. [12, 13] following the criteria that they must include integrated contributions which survive on large scales. 

Let us discuss this step in some details. From Refs. [12, 13] we are left to take the large-scale limit of the following 
expression 


VO 


(first — order)^ = / drj + iip'ip + — Ii{r]^)). 


(3.1) 


Jv* 

In Eq. (3.1) ri* is the linear fractional intrinsic temperature fluctuation on the last scattering surface, and /i(? 7 *) is 
minus the linear ISW effect, with 


/i(77)=2 



(3.2) 


while and are the first-order perturbation of the photon wavevector and background geodesic respec¬ 

tively [12, 13] 


fc(i)0(r;) 


-2ip + h{ri), 
rv 

2 / dfi[-(p + {ri- fi)(p '], 


-2 f drj [ipe^ + (?7 - , 


(3.3) 

(3.4) 

(3.5) 


where e* = —n* is the unit vector specifying the line-of-sight direction and we have dropped monopole terms due to 
the observer. In Eqs. (3.1)-(3.5) the integrands must be evaluated along the background geodesic (77, (ryo — ?7)e*), e.g. 
in Eq. (3.2) ip' = Lp'{fj,x. = —n(? 7 o — ?))), while a * denotes quantities evaluated at the last scattering surface, e.g. 
= 73(77*, X = -n{r]o - 77*)). 

Next, following Ref. [15] we notice that Eq. (3.1) can be further simplified. We use the relation 


and the integration by parts 


a.(i)o3,(i)*g, ^ _2 


2/ dr773'Ji(?7) =-[/i(r7*)]^ 

Jrio ^ 


(3.6) 


(3.7) 


Therefore Eq. (3.1) becomes 

rvo r 1 1 

(first - order)^ ^ + 73 *) - 473 * / d'q{'q^ - r])p, (3.8) 

7?* L 7?n J ^ -'77* 


^ The results in Refs. [13, 15, 42] have initial conditions corresponding to our a^i = 0. 



where we used the decomposition of the total derivative along the background geodesic for a generic function 
f{r],x^{r])), f = df/dr] = df/drj + dife^. In the following we will use the large-scale solution ri* = —2(p*/3. 

We are now able to provide an expression for the second-order ISW effect from Eq. (1.4). It is separated into two 
parts, the early ISW effect due to a non-negligible radiation component at last-scattering, and the late ISW due to 
the expansion growth suprresion in a A dominated epoch. Correspondingly we will split each integral into two parts 
as drj = drj drj, where rjm represents the epoch when full matter domination is reached. Before that 

time the early ISW effect is in action. Therefore we write the second-order ISW effect from Eq. (1.4) as 


with 


ATa 

T 


(n) 


ATa ATa 

-jr + -jr 


(3.9) 


1 ATa®”'^ 

2~Y~ 


dr] (4) -f 4')'(x, 7]) 


x=-n(r]o-v) 


(first - order)2^,iy, 


(3.10) 


where the (first — order)gg^j.jy corrections are obtained from Eq. (3.8) by keeping all those contributions which would 
vanish in the limit of full matter domination at the epoch of last scattering. We will discuss the early ISW effect in 
details in section IIIC. 

In the following we will focus on the late ISW effect 


1 ATa 

2^r~ 


di] (4> + 'I')'(x, ry) 


+ (first - order)f,,^tg 

x=-n{?7o-??) 


(3.11) 


using the results of Sec. II. Taking the time derivative of Eqs. (2.22) and (2.23) the late second-order Integrated 
Sachs- Wolfe effect is given by 


1 ATa 


2 T 


dr] 




2 ( -1 - -^(ani - 1) ) 9m9'{r]) + B[{r]) + 4g(?y)5'(?7) | ipl 


+ / di]g^B(r]) 


V ^ (SVoi^iTo) - 3V (^Voi^iTo) 


<.= -ii{rio-ri) 

-b (first - order)f,,tg .(3.12) 


x=-n(rjo-v) 


We recall that in order to obtain the expression (3.12) we have dropped from Eqs. (2.22) and (2.23) the terms 
[B 3 {r])\/~‘^did^{d^ipodj(po)-\-B 4 {r])d''(podig:>o], which give rise to the second-order Newtonian piece and give a negligible 
contribution on large scales, and in Eq. (3.12) the (first — order)j^^jg terms are obtained by taking that part of Eq. (3.8) 
which is due to late integrated effects 


(first - order) 



rv 

ifif' + / dfjif 






(3.13) 


For simplicity of notation in Eq. (3.12) we have introduced 


B(r]) 


4g'(ry)\ 
9m 3 gm J 


4 g'(rf)g(??) 

3 9l 



2 g^(g) 

3 g^ 


e'(g) ■ 


(3.14) 


Eq. (3.12) is the starting point to evaluate the CMB bispectrum induced by the late Integrated Sachs-Wolfe effect which 
arises from the explicit time-dependence of the linear gravitational potentials during the late accelerated phase [19]. 
A simple but important comment about Eq. (3.12) is that the whole effect is vanishing in the case of a vanishing dark- 
energy component. Moreover we can identify two fundamental contributions. One, proportional to [—5(a„i — l)/3], is 
directly sensitive to the primordial non-Gaussianity set in the early Universe from inflation. The importance of this 
term comes from the observation that in the case of strong non-Gaussianity (|ani| ^ 1) the ISW on large scales can 
be strongly amplified. The remaining terms are due to the non-linear evolution of the gravitational potentials after 
inflation and to second order-corrections in the temperature anisotropies. In the case of a high level of primordial 
non-Gaussianity one would expect these terms not to give the dominant contribution. However accounting for them 
is crucial in the computation of the second-order radiation transfer function for GMB temperature anisotropies, as 
we will discuss later. 
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B. Angular decomposition 

In order to compute the angular CMB bispectrum from the second-order ISW effect let us start from the expression 
of the (today) observed CMB anisotropies expanded into spherical harmonics as it is standard in the literature, with 
multipoles aim given by 


a,m = I d^n^^Y.Un). (3.15) 

In terms of the multipoles aim the linear and non-linear parts of the temperature fluctuations correspond to a linear 
Gaussian part and a non-Gaussian contribution 

a£m = a^m + a^m- ( 3 - 16 ) 

Before compnting for the second-order ISW effect let us recall briefly that at linear order the multipoles af^ are 
expressed as (see, e.g [43]) 


aL=47r(-I)'y (3.17) 

where (k) is the linear radiation transfer function which describes the relations between the initial fluctuations 
^ii(k) and the observed temperature anisotropies. In the following we take this initial value at the epoch of matter 
dominance ipuiY) = (pmik). For the linear Sachs-Wolfe effect acting on large scales (small k) one finds A^^^(fc) = 

ji{k{rio — ?7*))/3 where ji{x) are spherical Bessel functions of order £. For smaller scales one can compute 
with numerical codes such as CMBFAST [44]. Expressing the non-linear part in terms of the initial potential 
fluctuations corresponds to obtain a second-order radiation transfer function. The computation of the radiation 
transfer function at second-order in perturbation theory is very useful in that it allows, for example, to see what 
are the effects of the non-linearities (set on super-horizon scales) at the recombination epoch on different scales of 
the CMB anisotropy pattern. In this paper we are able to compute the second-order radiation transfer function 
corresponding to the large angular scales of CMB anisotropies. For small scales (big k) the computation requires 
to solve the non-linear Boltzmann equations for the photon-matter fluid. A comment is in order here. Notice that 
the primordial non-Gaussianity contribution proportional to [—5(ani — l)/3] will be transferred linearly to the CMB 
anisotropies, in that it is already a second-order contribution. Thus the crucial information to compute the second- 
order transfer function is carried by the remaining contributions in Eq. (3.12), which so far have not been properly 
taken into account.^ 

Let us sketch how to obtain afm for the late ISW effect (3.12) taking as an example just the term proportional to 
the initial non-Gaussianity A(n) = — dy ^(ani — 1) gmg'(v)‘PoU=-n(rio-T)) ^ ^llq. (3.12). Let us Fourier expand 
A(n) 

dryAW(ri)[</^^](k)e-*'^-fo°-fo, (3.18) 

(27r) 

where we set A = —10(ani — l)/3 and [i^Q](k) denotes the Fourier transform of If v^e now make use of the 

Legendre expansion + l)i^ jc{kx)Pe{k ■ x) we find 


JA r^o 

/ (^d^9^9'{v)je{Hvo-ri))[‘Pl]0^)Peik-n). (3.19) 

J r}m 


Thus inserting Eq. (3.19) into Eq. (3.15) we get the corresponding 


NL 


a,^^(A)=4^(-*)^ 


d^k 




rvo 


/ drj 


9'iv) . 


je{k{r]o - ry)) 


I"/m(k) 


(3.20) 


^ The analyses of the CMB bispectrum performed so far, as for example in Refs. [36, 37], adopt just the linear radiation transfer function, 
unless the bispectrum originated by specific secondary effects, such as Rees-Sciama or Sunyaev-Zel’dovich effects, is considered. 



10 


where we have used that ipm = Om^Po, the relation Y^’^(k)Y^m(h) = (2£+l)P^(k-h) and the orthonormality 

of the spherical harmonics. A [(^^](k) is a convolution 


J d^kid^k2S^^\'ki+k2-'k) 


-y(a„i - 1) 


p,n{kl)pra{k2) ■ 


(3.21) 


Following similar steps one can compute the contribution to from the other terms in Eq. (3.12). The extension 
of our example to most of them is straightforward, while for some terms in the (first — order)^^^^^ contribution (3.13) 
some care must be taken. The details of the computation can be found in Appendix B, and here we report the final 
result, namely 


NL 


= 47r(—z)^ 


d^k 

(^ 


(4^)^ E E 

LlMx L2^2 


'Ko{k)Af\k) + K,{k)Af\k) + K2{k)Af\k)] r/„(k) 
d^ki d?h 


mMi M 2 
L 2 


(27r)3 (27r)3 


Pmikl)pm{k2) ALjL2(fci,fc2) Tl^Mi (kl)yL2M2 (k2) , 


(3.22) 


where Knik) are convolutions in Fourier space expressed in terms of some kernels /„(ki, k 2 , k) as 

1 


Kn{k) = 


(27r)3 


d^kid^k25^^\ki + k2 - k)/„(ki, ka, k)v3m(ki)v3m(k2 


with 


/o(ki,k2,k) = --(oni-l)-!, /i(ki,k2,k) = 1 


/2(ki,k2,k) = 3 

k given by k = ki + k 2 , and correspondingly 


(ki-k)(k2-k) ki • k2 


k^ 


fc2 


= 2 r -■ 

J r?. 

^rik) = f 

J n 


dr]^-^j(,{k{r)Q - 77)), 

r/m 9m 


Vo T^l („\ 

dv^Mikim - ri)), 

9m 


Af^\k) = - dTjB{r])jj,{k{Tjo - Tj)). 


In Eq. (3.22) = f d^nY£m(n)VLiMi (^M 2 M 2 (n) is the Gaunt integral and 


^LiL2{kl,k2) = -4 
+ 2 


I dr] 

Vm 9' 


9"{V) . 


-9iv) . 


jL^ikl{r]o - r])) / dfj - jL2ik2ivo - v)) 

Jtjo 9m 


[ dr] jhAkiim v)) 

Q'iv) 1 

2 / dr] jL 2 {k 2 {r]o v)) + JL 2 ik 2 {vo V*)) 

'rim 9 m 

Jrim 9m 0 


(3.23) 

(3.24) 

(3.25) 

(3.26) 

(3.27) 

(3.28) 


(3.29) 

Eq. (3.22) is clearly the generalization of the linear relation (3.17) with the functions A^^‘^\k) and Ai^^L^{ki, k 2 ) 
playing the role of coefficients of the second-order radiation transfer function relating quadratic curvature perturbations 
to the observed temperature anisotropies. ® For the second-order ISW effect the transfer function depends on the 


^ In fact for the last term of Eq. (3.29), Jl 2 (^ 2 ('^?o — '^*))/3, corresponding to the last term of Eq. (3.13), we are assuming full matter 
domination at last scattering 77 *, so that we used (/?* (/? 7 n- However in the final formulae we display in Sec. V we account also for the 

small differrence between (/?* and (prn- 

® Notice that it is possible to rewrite Eq. (3.22) with a more general formula as 


<^z = (4-)" E E 

L\M\ L 2 M 2 


d^ki (fik2 

(27r)^ (27r)^ 


(kl )v^ m (k2) 


^LiL2^^~^A2) 'll,iMi(kl)yL2M2(k2) , 


(3.30) 
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growth suppression factor g(ji). As we anticipated, for the initial non-Gaussianity, parametrized by Oni and contained 
in the quadratic term KoCk), the transfer function is exactly the same as in linear theory, while the transfer 

function becomes non trivial for the remaining terms quadratic in the potential perturbations. In particular the part 
of expressed in terms of the function AL-^^L^{ki,k 2 ) comes from the last three terms in Eq. (3.13). Notice that 
for these contributions the dependence of AT/T on h does not enter only through a single angle k • n, as it is usual 
in linear theory. 


C. Calculation of the early Integrated Sachs-Wolfe effect at second order 

At the epoch 77 * of last scattering the universe is still not completely matter dominated, and the residual radiation 
component makes the potentials decay in time. Such an effect gives a non negligible contribution to the large-scale 
CMB anisotropies, since it almost mimics a Sachs-Wolfe effect. For example at the linear level (for a vanishing A) 
the early ISW effect gives a correction to the power spectrum normalization of the order of 20% [26]. A simple 
estimate from Eq. (3.8) shows that also at second-order the early integrated effects must be taken into account giving 
a significant correction to the second-order Sachs-Wolfe effect (1.2). 

We consider a universe filled by radiation with energy density and pressureless matter with energy density pm, 
supposing that the cosmological constant term is negligible at the epoch of last scattering. In this case the linear 
growing mode of the gravitational potential is 


with [53, 54, 55, 56] 


where 


ip = 


7^* 


Fiji) 


F{p) 


H r ^ _ 3 


2 8 16 
1% ~ IV ~ IV 


16^1 -I- y 
15y^ 


Pm 

y= — oc a. 
Pi 


(3.32) 


(3.33) 


(3.34) 


H{a) is the Hubble expansion rate, 07 is some epoch when radiation is dominating over matter and the * indicates 
evaluation at the time of last scattering. Let us briefly recall how to compute the linear early Integrated Sachs-Wolfe 


effect 


= 47r(-i)' 


(Pk 

(^ 


dp'p'ik, T])je{k{r]o - rj)) 


YLik), 


(3.35) 


where the time derivative of the gravitational potential is due to the transition from a residual radiation component 
to full matter domination reached at the epoch pm- Since most of the contribution to the early ISW effect comes from 
near recombination at 77 * the integral (3.35) can be computed by evaluating the Bessel function at leaving [26] 


rrim 

2 J dp(p'{k,p)f=i2(p{k,p)\^’^je{k{po-p^))=2<p^(k)—je{k{po-p^)), 

where in the last line we used Eq. (3.32) with AF = F{p)\'^’^. 

Let us now consider the effect at second-order. The early ISW effect reads 


AT 2 

~Y' 


rVm 


+ 'Sy{x,p) 


+ (first - order)2^,i 

x=-n(?7o-??) 


(3.36) 


(3.37) 


with 


Al,i,(kl,k2) = 


(_l)£+ii+l ,2 1 y /„(ki, k 2 , J drr^jLiikir)jL^{k2r)j({kr) 


-I- (kiTi) ■ 


(3.31) 
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where the (first — order)gg^j.jy corrections are obtained from Eq. (3.8) by keeping all those contributions which would 
vanish in the limit of full matter domination at the epoch of last scattering 


drj 


ipLp + Lp 


(first - order) 2 ^^iy = -4 / 

rVo 

- 4(/?(. / ^ 77 ( 77 * - r])ip, 


dfjip 


+ 2^U'nm]ri*) +/i( 77 ,„; 77 *)/i( 77 *) - -<p*Iiir]m]V*) 


(3.38) 


where we have splitted defined in Eq. (3.2) as 

rm I'Vm rm 

hiv*) = -‘^ dr](p'=-2 dr](p'-2 drjip'= Ii{r]m;r]*) + hiVm) ■ 

J n* d rim. 


(3.39) 


Notice that, for example, the last term in Eq. (3.38) proportional to would give a vanishing contribution in the 
limit of matter domination when the gravitational potential Lp is constant in time. 

For the moment let us focus on the first term on the R.H.S of Eq. (3.37). The corresponding multipoles are easily 
computed as in Eq. (3.35) 




d^k 

(^ 


rVm 


dT]{<^ + d')'(k,77) jiikirjo - 77 )) 




(3.40) 


with the standard approximation 

rVm 


1 r'^ 1 

- J dr]{<^> + d')'(k, 77 ) je{k{r]o - 77 )) ~ -($ + jeikim - V*)) ■ 


(3.41) 


Therefore what we need is the expression for the second-order gravitational potentials around the epoch of last¬ 
scattering when both radiation and matter are considered. The details of such a computation are reported in Appendix 
C, where we consider the evolution of the second-order gravitational potentials for a universe filled with radiation and 
matter. 

First of all we find a relation between the gravitational potential $ and 4' 

^ = '^ + Alp^-Q, (3.42) 

where Q can be written as 

2 = (3-43) 


with K. given by Eq. (1.3) and 


g = -4 


1 + y 

4-k3y 


+ 2FF - 


16 + 5y 
2 l + y 



(3.44) 


In Eq. (3.44) y = Pm/P-y is the ratio of the matter to the radiation component, and a dot stands by a logarithmic 
derivative with respect to the scale factor, ' = d/d In a. Eq. (3.42) is the analogue of the relation found in in Ref. [25] 
for a matter-dominated Universe and of Eq. (2.11) (or Eq. (A.11)) for the case of matter plus a cosmological constant. 
The evolution of the gravitational potentials is then derived from the equation for 4' 




where AC 2 is given by Eq. (2.14) evaluated for a system of radiation plus pressureless matter 


® T 

-47 

' =+ 

+2 + q-Ea(2 

VVp \ 


2 p 


EaC 2 = 2- \2RF‘^ - 2(2 - R)F2 - 4(1 - R)FF 
‘^P Fi V 


vl, 


and for simplicity we have introduced the function 


(3.45) 


(3.46) 


R{y) 


1 + d 
4-k37/ 



) . 

4 + 37/; 


(3.47) 
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Eq. (3.45) has been obtained by using the previous relation (3.42) in the (O-O)-Einstein equation on large scales 
combined with the expression for the second-order curvature perturbation (2.13) (see Appendix C). 

In the case of adiabatic perturbations, for which C 2 = constant, the integration of Eq. (3.45) yields 


a Ja, Vp 




Q-i-AC2 

2p 


C^. 


(3.48) 


We thus find the usual linear relation plus a source term quadratic in the first-order perturbations. The last term is 
a decaying mode that can be neglected. 

Finally, using Eq. (3.42), we obtain 


= -2F{r]) C 2 + 
^ 'Q-2 


10E2 


1 

V~p 


+ 2 ^ J ((1 - 2i?)E2 -h 2(2 - R)F^ + 4(1 - R)FF^ 




da 

7^p 


Q 


JC, 


(3.49) 


which is the expression to be used in Eqs. (3.40) and (3.41) when eavluating the difference (T -|-<I))|^™. The conserved 
quantity C 2 is determined by the primordial level of non-Gaussianity parametrized as C 2 = 2aniCi 7 with Ci = —ip^jF^. 

To compute the multipoles for the additional terms (3.38) one has to follow the same procedure described in 
Appendix B. We will use Eq. (3.32) for the epoch when radiation is not negligible (i.e until the time r\ra) and we will 
adopt the growth suppression factor g[r\) as in Eq. (2.3) from full matter domination onwards, with 


9m 5 

-T * 


(3.50) 


where denotes the value F{rjm)- Moreover we will use the evaluation (3.36) for the linear early ISW effect. 

Thus we are able to calculate the radiation transfer function for the early ISW effect from Eq. (3.49) and Eq. (3.38). 
We find 

d^k 


az = 


(27r)^ 


(4^)^ E E 


Xo(k)Af ^ft) + + X 2 (k)Af \fc)] y/„(k) 

d?ki d?k2 


Ml M 2 


LiMi L 2 M 2 


(27r)3 (27r)3 


(p*(ki)v3H<(k2) ALiL2iki,k2) TLjMi(ki)yL2M2(k2) 


(3.51) 


where (/?*(k) is the gravitational potential at the epoch of last scattering, the convolutions Ar„(k) are given by 

1 


i^n(k) = 


(2^)3 


d^kid^k 26 ^^yki + k2 -k)fn{ki,k2,k)(p^(ki)p^{k2), 


(3.52) 


with the kernels /„(ki,k 2 ,k) defined as in Eqs (3.24) and (3.25), and the transfer functions are given by 


Af^\k) = ji{k{r]o - V*)) 


. 1 ( 2 )/ 


(3.53) 


1 




— J — ((1 - 2R)F‘^ -k 2(2 - R)F‘^ + 4(1 - R)FF^ 


jiikirjo - 77*)) 


4 AE ,, , 

+ Mkivo - V*)) 


^r\k) = ^ 


1 


Q Vp 

2 a . 


rE 


Q 


je{k{r]o - r]^)), 


(3.54) 

(3.55) 


and 


AL^L2{ki,k2) — ( 2 ^2 + 


AE2 


■ (Ul \\ 

■jL2[k{vo-W 


F, 


AF\ AF g'{r]) 

jLiikim - r]4) jL 2 ikir]o - V*)) + 4:—jL2{kirio - V*)) / dpc - jhAkipo - v)) 

J J rim 9m 

pm pp\ po gfp 

/ dr]{r]^-r])—^jLAkiilo-V*))+ dg{T]^ - rj) c -jLi (fc(?7o - ??)) 

Jr;. Jrim 9m 


F 4 




drj 


jL2ikivo - v)) 


'V* 




dfj 


Fig) 


Fjm 

m 


9{v) . 


P -jLr{k{go-fi)) + / dgc -jLi (^(r/o - ?))) 

^ * .1 'n 9m 


'rim 


(3.56) 
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where c = Fm /and 


AF = F{7i^) - FiTj,). (3.57) 

From Eq. (3.53)-(3.56) we see that the second-order early ISW effect can be approximated on large scales as a 
Sachs-Wolfe effect (whose transfer function is porportional to je{k{r]o — ry*))). 


IV. CMB ANISOTROPIES FROM SECOND-ORDER TENSOR PERTURBATIONS 


The second-order tensor contribution to the CMB temperature anisotropies is given by the integrated effect [13] 


AT2 

T 


(n) 


1 

2 



dvxif (x, 77 ) 71 * 71 -^ 


x=-h(r]o-v) 


(4.1) 


On large (superhorizon) scales the tensor modes remain constant (see, e.g., Ref. [45]), while they start to evolve when 
they are reentering the cosmological horizon. Therefore the main effect on the CMB anisotropies comes from late 
times, and therefore we can take 77* « ijm- Notice also that the second-order tensor modes are determined by a 
contribution produced during inflation and an additional contribution generated by the post-inflationary evolution. 
Since the level produced during inflation is very small (of the order of the slow-roll parameters) compared to the 
post-inflationary contribution, as shown for example in Ref. [24], we will neglect it in the following. 

For a nonvanishing cosmological constant A and pressureless matter we show in Appendix A how to obtain the 
evolution equation for tensor modes 


X2, 


(4.2) 


where the source term is given by 

f) = + 8 V -2 - d^dkR^j - d^djR\) + , (4.3) 

and the traceless tensor 


^ (V(/?)^ 6i + AnGa'^p (^vivik - 
2E^iz)p{n^)\ W, 


= 3^1 + 


aVo9fc(/?o - ^ (V(/?o)^^i) ■ 


3flom (1 + z) 

It proves convenient to rewrite the source term as 

f) = s{p)T){-k ), 

where 

, , o 2, ^ ^ , 2E^z)p{n^) 

= -8s (,) (_i + + > 

and 

V'^Tij = Sij = V'^&oSij + didjQo + 2 {didjpo'^^y^o - didk(pod’'dj(po) 


(4.4) 


(4.5) 


(4.6) 


(4.7) 


with 00 defined in Eq. (A.31). 

Let us first consider the solution of the second-order tensor modes. We can Fourier transform X 2 j 9-^*^ source 
term T* and decompose them into the ct = -I-, x polarization modes. 


/ fF'k 

(^e*‘""A2,(k,77), 

X2j(k,77) = x+(k, 77 )e+f(k)-bXx(k, 77 )e''}(k), 


(4.8) 

(4.9) 
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where (k) are the polarization tensors, and similarly for the source term. It is then easy to show that the solution 
of Eq. (4.2) can be written as 

X^,(k,r;)=r;.(k)Mfc,r7), (4.10) 

where 

h{k,'n) =xi{k,v) [ - X2{k,ri) f , (4.11) 

and xi X 2 the solutions for the linear gravity-waves amplitude 

x'' + 2Hx' + kX = 0, (4.12) 

= X 1 X 2 — X 1 X 2 being their corresponding Wronskian. In Eq. (4.10) ir®^ (k) represents the amplitude of the second- 
order tensor modes fixed by the source term (a quadratic combination of linear scalar modes). 

Let us now consider the expression for the multipole coefficients aim for tensor perturbations. We follow some 
standard References [46, 47, 49] where the reader can find more details on the calculations. From Eq. (3.15) we can 
define for a given k-mode the coeffiecients 

aUk)= I (4.13) 


where we have expanded in Fourier space the anisotropy AT{x,n,ri)/T = (27r^) AT(k,n,r/)/T to be 

evaluated today and at the origin. Thus from Eq. (4.1) we get for the 

'1 r 




(k) = (-1)^ 


d?7X+(k,?7) 




e+- 


polarization state, 

(k)n*n^ y] /(2£' l)je'ik{r]o - r]))Pe>{h) 


(4.14) 


where we have used the Legendre expansion = '^£{2£ + l)i^ je{kx)Pe{k ■ x), and a similar expression holds for 
the X polarization state. 

It now proves convenient to perform the angular integral for a specific k-mode, by choosing the coordinate system 
such that z = k. Statistical isotropy of observable quantities, like the angular power spectrum or the angle-averaged 
bispectrum (and higher-order correlation functions) ensures then that we can take the sum over the different modes in 
k space at the end. Alternatively we will also display a general formula for a generic k (and for the coefficients aim)- 
With such a coordinate system choice the integral in Eq. (4.14) can be evaluated exactly as shown in Refs. [46, 47, 48] 
to give 


aL(k) = 


X 


(^_iYy^yq2lTT)^J^^^{6m2 + Sm,-2) j drix!+{k,vi) 

je +2 _^ 2 _ il _ I_ di -2 

_{2i + 3)i2i+l) {2£ + 3)i2i-l) {2i + l){2i - 1)_ 


- a / {2£ -X 1 ) 


- I {£ + 2)1 

-T, -^ 


(Y 


-1- A 


\ n-' 71. \ 



dr] h'{k, rf) 


jejkivo - v)) 

[k{m - v)V 


(4.15) 


where the spherical Bessel functions are evaluted in k{rjQ — rj). In the last line we have used Eq. (4.10), 7+(k) being the 
-I- component of the amplitude Tj (k), and we have applied the recursion relation ji{x)lx = {ji-i{x) + ji+i{x))/2£+l. 
For the x polarization state one has to replace { 6 m 2 + dm,- 2 ) with i{Sm ,-2 — dm 2 )- 

Following Ref. [48], one can also give a formula holding for a generic k-mode, by using the Wigner-D functions [50] 
to perform a rotation of the spherical harmonics. One thus finds 


a-tm — {~i)^~^\/ { 2 £ + 1)^ 


'(7 + 2)! 


, (^-2)! J 

+ (x+(k,r7)+iXx(k,»7))-DL,-2(k) 


J rim 


[Km - vW L 


ix'+{Kv) - ix'x{Kv))D'm2{^) 


(4.16) 


where D\^^i (k) are the Wigner-D functions corresponding to a rotation i?(k) which brings k from (0,0,1) to a generic 
direction. 

The evolution of the second-order tensor perturbations is given by h{k, rj) having applied the Green method. The 
solutions for the corresponding homogenous equation (4.12) (which is the one for linear tensor modes) are well known 
for the limiting cases flm 1 and Da —> 1, and correspondingly analytical expression exist for xhj (see, e.g., Ref. [28] 
and [51]). However the full solution requires a numerical evalutation of Eq. (4.12). 



V. SECOND-ORDER RADIATION TRANSFER FUNCTION FOR LARGE-SCALE CMB 
ANISOTROPIES: SUMMARY TABLE OF RESULTS AND CONCLUSIONS 
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In this section we finalize our computations by providing a summary of the formulae for the full radiation transfer 
function for the CMB anisotropies on large angular scales due to the second-order Sachs-Wolfe effect in Eq. (1.2), the 
(early and late) second-order ISW contribution (3.12) and the second-order tensor modes, Eq. (4.1). The section can 
be regarded as a summary of the main results of this paper. 


A. Large-scale CMB anisotropies from second-order scalar perturbations 


The second-order large-scale CMB anisotropies can be written as 


aZ = 


dPk 

(27r)3 


(4^)' E E 

L 2^2 


Koik)Af\k) + K^ik)Al^^\k) + K2ik)Af\k)] Y/^(k) 
'■ d^ki cPk2 


Ml M 2 

L 2 


(27r)^ (27r)3 


(p*(ki)v3H.(k2) AL^L2iki,k2) TLjMi(ki)yL2M2(k2) 


(5.1) 


where (/?*(k) is the gravitational potential at the epoch of last scattering, A„(k) are convolutions in Eourier space 
expressed in terms of some kernels /(ki, k 2 , k) 


Kn{k) 


1 

(^ 


(fkid^k2S^^\ki + k2- k)/„(ki,k2,k)(/?*(ki)v3H.(k2), 


(5.2) 


with 


/o(ki,k 2 ,k) 

/i(ki,k2,k) 

/2(kl,k2,k) 


5 

~ ^(Onl ~ 1) ~ 1 ) 

1 , 

(ki ■ k)(k2 ■ k) _ ki ■ k2 

k* k"^ 


(5.3) 

(5.4) 

(5.5) 


k given by k = ki -|- k 2 . 

In Eq. (5.3) Oni is a parameter which specifies the level of primordial non-Gaussianity generated in the early Uni¬ 
verse from inflation (or some other alternative mechanisms for the generation of the cosmological perturbations) [8]. 
The functions A^^‘^\k) and AL^L^{ki, ^ 2 ) represent the second-order radiation transfer functions, and they are given 
by the sum of the transfer functions for the Sachs-Wolfe effect (1.2) and the Integrated Sachs-Wolfe effect (3.12). 


• Second-order Sachs-Wolfe effect 


Let us first consider the Sachs-Wolfe effect. Notice that the Fourier transform of 3/C/lO corresponds to 
the convolution A' 2 (k) defined by Eq. (5.2) and Eq. (5.5). Thus it is straightforward to compute the corre¬ 
sponding multipoles for the Sachs-Wolfe effect by using the same steps described in Section III. From Eq. (3.15) 
we find 


Af^\k) = ^je{k{r]o-V*)), (5.6) 

= ^j£{k{r)o-r]^)), (5.7) 

Af'^\k) = -ijf(fc(r?o - ?7*)), (5.8) 

AL.L^ikiM) = 0. (5.9) 
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• Second-order Early ISW effect 

In Section III C we found 

jeiHm - V*)), 


am 

^ yl £ ^(^(1-2R)F^ + 2{2-R)F^+4{1 -R)Ff) - V*)) 


Af Hfc) = I - V ’ 


(5.10) 


(5.11) 

(5.12) 


AF2 2 AF 


3 E* 

■ rm 

, g'ig) 

dr]c - 

J V-m 

gm 

r'nm 

, F"{g) . 


jhAHm - v*))jLAk{m - v*)) 

1 AF 

(Hvo - v)) -E^jL2 {kim - V* 

r ,(c 

r fv^ p(fj) 


r° giv) 

/ drjc -j 

fr^m 9m 


f dr){r)^-r))^^jLAk{vo-V*)) ^ jL^ikivo - V*)) 

Jr}^ * J * 


E* 

F'iv*) 


( d?7(??H.-??)c—jLi(fc(?7o-??)) ^-^^JhAkivo-g*))- 

J rim gm \ 


(5.13) 


• Second-order Late ISW effect 

From the results of Sections 11, III A and IIIB we have 


A°/‘^\k) = 2 / d-qc^^-^jiikirto-r])), 
Jrim g^ 

Al^‘^\k) = j dr]c^^dpljj,(^k{r]o-r])), 

J r}m 9m 

Af‘^\k) = - dr]C^ B{r])je{k{rio-v)), 

9r)m 


(5.14) 

(5.15) 

(5.16) 


r° Q”iv) r q(v) 

ApL 2 {ki,k 2 ) = -4 / dr] c-—- jLi{ki{r]o - g)) df) c (fc 2 (r?o - ??)) 

Jrim 9m JriQ 9m 

Q^iv) p'iv) 1 

+ 2 dr]c^AJlj^^(^ki{r]o-g)) 2 dr] c^—A j - g)) p j L^{k^{gfi 

J n-rrr gm 

J Tlrrr g-m o 


g*) ■ 
(5.17) 


Let us recall the definition of the different symbols and quantities appearing in the previous formulae. 

- Symbols 


g 


conformal time 
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a scale factor 

Qi epoch of radiation domination 

rj^ epoch of last scattering 

rjm epoch of full matter domination 

ryo present epoch 

(/?* (linear) gravitational potential at the epoch of last scattering (5.18) 

• = ^ 
din a 

' = - 
dr] 


- Functions 

• F{t]) gives the linear growing-mode solution for the gravitational potential (p for adiabatic perturbations of 
pressurelss matter and radiation fluids [53, 54, 55, 56] 


where 


with 


F{r]) 


3 2 8 

5 ~ Tv 


16 mVTTy 
I5y^ 15y3 ’ 


y = — (X a . 
Pi 


(5.19) 


(5.20) 


c = Fm/F^ 


(5.21) 


and 


AF = F{r]m) - F{ri^). 


(5.22) 


• g(ri) is the growth-suppression factor for a flat Universe filled with a cosmological constant A and a pressureless 
fluid whose exact form can be found in Refs. [29, 30, 31]. In the A = 0 case g = 1- An excellent approximation 
for g as a function of redshift 2 : is given in Refs. ]29, 30] 


goi^m 


^rrJ ~ + (1 + flm/2) (1 + Ua/70) , 


(5.23) 


with Urn = + z)^/E‘^{z), Oa = ^oa/E‘^{z), E{z) = {I + z)Ti.{z)/H q = [r2om(l + ^)^ + J^oa] and 

riom, floA = 1 — flom, the present-day density parameters of non-relativistic matter and cosmological constant, 
respectively. We have normalized (/(ry) so that at the present epoch g{z = 0) = 1. The quantity g^ corresponds 
to the value of g{r]) at the epoch rjm when full matter domination starts. A good approximation is 


9m ~ 



(5.24) 


• The function Q enters in the difference between the two second-order gravitational potentials 4' and $, Eq. (3.42), 
and it is given by 


Q = -4 


1 + d 

4 + 3j/ 


• 9 • 1 6 -l- 5?/ 

F^ + 2FF + -i 

2 1 + y 



(5.25) 


• The function R{y) has been introduced for simplicty and it is defined as 


R{y) 


i + y 
4 + 3yy 



) • 

4 + 32/7 


(5.26) 
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• The functions Bi{r]) and i? 2 ('' 7 ) enter in the large-scale solution of the gravitational potentials T and 
Eqs. (2.22) and (2.23), 


where 


Biiv) = Ho ^ (/o + 30o™/2) ^ 

rjm 


r'n 


^ 2 ( 77 ) = 2Ho ^ (/o + 3flom/2) W dr)H^(? 7 )| 2 (/(? 7 ) - 1)^-3-f 3 fl™(? 7 )|C'(??, 7 )), 

' Vm 




C{v,fi) = 9 (j7)a(?))|5(77)H(?))-g(?))^7;;^H(?7) 


\9) 




6 ( 77 ) = 


Piv) 


Bijj) is defined by 


BP) = 


( 77 ) 

^2(77) 4 5'(?7)\ 4g'(?7)5(7?) 


3 5r7 


3 9?, 


3\ 25 ^( 77 ) , 

Pv) + 0 + o^—e (77 ■ 

2/ 3 gi 


(5.27) 

(5.28) 

(5.29) 

(5.30) 

(5.31) 

(5.32) 


B. CMB anisotropies from second-order tensor modes 

In Section IV we found that the contribution of a single k-perturbation mode, with the coordinate system choosen 
with z = k, is given by 


i(k) = i-iY^VWTT) 


+Sm- 2 ) B+ik) j dgh\k,r]) 
(k) = i-iY^PiWTT)^ zi 6 m ,-2 - Tx (k) dg h'{k, 77 ) 


Jt(fc(7?0 - 77 )) 
[fc(770 - 77 )] 2 

jeikjm - 77 )) 

[M7?0-7?)]2 ’ 


(5.33) 

(5.34) 


for the + and x polarization states respectively. 


71|_(k) and irx(k) are the -I- and x component of the amplitude of the second-order tensor modes, H* (k) = 
/ d^y e-^y '^T){y) where 


with 


+ didjOo + 2 {didj(poV'^(po - didk(pod^djy^o) , 


V"0o = - 2 ((VVo)' - dA^od^d'^ipo) . 


h{k, 77 ) describes the evolution of the second-order tensor perturbations 

Hk,v) =Xiik,r]) [ dfi^P%^s{fj) -X2ik,g) [ dfj^P%^s{f)), 


W 


W 


where 


s( 7 ?) = -85^(77) 1 + 


2 E^(z)p(n^) 
3fIom(l + 

and xi and X 2 are the solutions for the linear tensor modes 

x" + 2nx' + k^x = 0, 

^ = X 1 X 2 ~ X 1 X 2 being their corresponding Wronskian. 


(5.35) 

(5.36) 

(5.37) 

(5.38) 


(5.39) 
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C. Conclusions 

In this paper we have computed the full second-order radiation transfer function for Cosmic Microwave Background 
anisotropies on large angular scales in a flat universe filled with matter and cosmological constant A. Our calcu¬ 
lation includes the second-order generalization of the Sachs-Wolfe effect as well as of (both the early and the late) 
Integrated Sachs-Wolfe effects, and is valid for a generic set of initial conditions specifying the level of primordial 
non-Gaussianity. We also accounted for the contribution of second-order tensor modes. Our results deal with large 
angular scales and represent the first step towards the computation of the full second-order radiation transfer function 
for CMB anisotropies. Its computation on small angular scales involves several non-linear effects, such as, for example, 
gravitational leasing, Shapiro time-delay, Rees-Sciama effects, and of course, dealing with the non-linear Boltzmann 
equation for the photon-matter fluid. Such computation will be the subject of a future publication [19]. 
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APPENDIX A: PERTURBATIONS OF A FLAT ACDM UNIVERSE IN THE POISSON GAUGE. 

1. Background equations 

Let us consider a flat Friedmann-Robertson-Walker (FRW) universe filled with a pressurelss fluid of energy mo¬ 
mentum tensor = pu^Ui, plus a cosmological constant A. The Friedmann equations are given by 

= ^(SttGp + A), (A.l) 

p' = -37fp, (A.2) 

where a{r]) is the scale factor, a bar denotes the background quantities and Ti. = a'/a, with a prime standing for a 
derivative with respect to conformal time rj. Another useful relation is 

2H' = -n^ +a^A. (A.3) 

Finally let us recall the definitions of the density parameters flm = SttGq^ p/{3TC^) and Oa = a^A/(37f)^ for the 
non-relativistic matter component and the cosmological constant respectively, with 12^ + = 1- 


2. Linear perturbations 

We briefly report the results for the linear perturbations (for a detailed analysis, see e.g. Ref. [57]). The perturbed 
metric tensor is defined by Eq. (2.1) and we will adopt the Poisson gauge for which LOi is a pure vector, i.e. diU)^ = 0, 
while Xij is a tensor mode {i.e. divergence-free and traceless d^Xij = Oj = 0 )- The matter component has mass 
density p{x.,ri) = p(? 7 )(l -I- d) and four-velocity = (Sq + v^)/a, with = —1. Each perturbation is expanded 

into a first-order (linear) and second-order part, e.g. for the mass density perturbation 6 = 5i(x, rj) + 52(x, ri)/2. The 
traceless part of the {i — j) Einstein’s equations in the Poisson gauge imply that (pi = 'ipi = ip, while its trace gives 
the evolution equation 


p" + 3Hip' + {2H' + n^)ip = 0 , (A.4) 

which, with Eq. (A.3), brings to Eq. (2.2). Vector modes are not generated from standard mechanisms for cosmological 
perturbation, such as inflation. Tensor modes obey an evolution equation obtained from the traceless part of the {i- 
j)-component of Einstein equations 


xl/' + 27fxl/-VVi, =0. 


(A.5) 
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However, for the purposes of the present paper, we can neglect the linear tensor modes in the following. 

The energy and momentum constraints provide the density and velocity fluctuations in terms of Lp (see, for example. 
Ref. [25] and [15, 28] for the A case) 


3. Second-order perturbations 


a. Scalar perturbations 


Let us first consider the second-order scalar perturbations of the metric in the Poisson gauge <j >2 = ^ and 11^2 = 

We will now show how to obtain the evolution equation (2.7) for the gravitational potential and the relation (2.11) 
which represent our master equations. 

The second-order perturbations of the Einstein tensor can be found for any gauge in Appendix A of Refs. [8, 23]. 
The perturbations of the energy-momentum tensor up to second order in the Poisson gauge have been computed in 
Ref. [25] for a general perfect fluid (see Ref. [8] for expressions in any gauge). For the pressureless case one has 

T 0 = — p ^1 -I- i5i -I- -62 + , 

T*o = —p v\ +-V 2 + {if + 5i) v\ , 

= p Vii + —V 2 i + —U} 2 i + { — ^(p + 61 ) Vii , 


T) = p , (A.8) 

where vf = v\vii. Note that the second-order velocity V 2 i is the sum of an irrotational component which is the 
gradient of a scalar, and a rotational vector which has zero divergence, = 0. 

We thus obtain the (f-j)-component of Einstein equations 

T" -b H{2'S' + $') -b ^V2($ - ^') -b {2H' + H^)<P - 4(77' -b H^) - SHpp' - 3 (d^pd^ip) - 4(pVV <5) 

- T) - ^nid^uj2j + djujl) - i(5*W2i' + djUJ2) + \ xk” + + 2d^pdjp + Apd^djp = 

SnGa^pvlvij . (A.9) 

We now proceed as follows. We take the traceless part of Eq. (A.9), S 2 G^j — \ 52 G^f. 5 '‘j = 87rG(52Rj — ^^ 2 T\ 6 j), 
namely 

- (T - $)-b ^ (V(p)^-b S^j + ld^dj{'S-<l>) + 2d^pdjP + 4pd^djp 

-i (5 *u;2,/ + dM') - i (x^/' + 2Hx^; - v2x* 2,) = S^Ga^p (A.IO) 

and we apply the operator did^ in order to get rid of second-order vector and tensor perturbation modes and to solve 
for the combination (4> — $). We find 


where Q is defined by 


T - $ = -4(p2 -b Q , 


\/2q= -P + 3N, 


(A.ll) 


(A.12) 


P = 

^ ^ I 1 


(A.13) 
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where 


and the quantity N is given by 


P*- = 2d^ipdjip + SirGa^pvlvij , 


V^iV = . 


(A. 14) 


(A.15) 


The combination — $) will be used in the trace of the (i — j) Einstein equation that we write in the form 
4'" + + a^A4' = Hi'S' - $') + a^A(4' - $) + - 4>) + da^At^^ + + hd^pd.ip) + SHpip' + 


'-22 
-pa , 


(A.16) 


where we have used Eq. (A.3). Inserting Eq. (A.11) into Eq. (A.16) and using Eq. (A.12) and (A.13) we find 

«'" + 3^4-'+ a^Ad- = HQ'+ a‘^AQ + N + p'‘^ (A.17) 
The next step is to write explicitly HQ' -\- a^KQ and N through Eqs. (A.12)-(A.15). Notice that by using Eq. (A.7) 


one can write as 


- d'p'djp' + H{d''pdjp)' + — a?h)d'pdjp 

P V ^ 


n- = 4 ^^ yd'p'd.p' + H{d\d,p)' + -(hW - a^k) 
where we have also used the background equation (A.l). Thus from Eq. (A.15) 


= 


AirGa'^p 


d\'djp' + H{d^pdjp)' + -{5H^ - a^A)d^pdjp 


(A.18) 


(A.19) 


Eor the combination HQ' + a'^AQ, with the definitions (A.12) and (A.15), we can write 

V^V^iHQ' + a'^AQ) = -W‘^{HP' + a^AP) + 3d^d^{HP'' + a^APM . 


(A.20) 


Using the background equations (A.l), (A.2) and the evolution equation (2.2) for the linear potential p we get 


PP = - ^ ^ [iHd^p'djp' + a^AHd'^pdjp + a?A{d'^p'djp + d'^pdjp')] , 


(A.21) 


and + a^AP*^ = 2a^Ad^pdjp — Ad^'p'djp'. We then find 

V‘^V‘^{HQ' + a^AQ) = 4 {dip'd'p') - 3did^{d^p'djp')] - 2a^A [^‘^{d.pd^p) - 3did^{d^pdjp)] . (A.22) 

We can use the relation f{flm) = d\nD+/d\na = 1 + g'{v)/{Hg{p)) and Eq. (2.3) to write 

HQ' + a^AQ = 2g^H^nm (2 ^^~ _ A + 3 ^ (^-^(d'pod.po) - 3V-'^d^d,id,pod^po)] , (A.23) 


and 

N = ^g^ + 1) V-^d.d^d^pod.po ), (A.24) 

where V“^ is the inverse of the Laplacian operator. Therefore, inserting Eqs. (A.23) and (A.24) into Eq. (A.17), we 
finally obtain Eq. (2.7). 

Notice that, in a similar way to N, we can write the quantity P defined in Eq. (A. 13) as 

p ^ 4^2 ^m) 0 d^pod^po ■ (A.25) 

Thus applying the operator to Eq. (A.11) and combining Eq. (A.25) with Eq. (A.24) we find the relation 

between 4* and 4) of Eq. (2.11). 
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b. Vector perturbations 


Let us consider here the second-order calculation which leads to the generation of second-order vector and tensor 
modes. We will strictly follow the findings of Ref. [28]. For the A = 0 case, this problem was originally solved in Ref. 
[42] starting from the results of second-order calculations in the synchronous gauge [58, 59] and transforming them to 
the Poisson gauge, by means of second-order gauge transformations [60]. In Ref. [28] the problem is solved directly in 
the Poisson gauge accounting for a non-vanishing A term. For a non-zero-A also Ref. [15] studies second-order vector 
and tensor modes in the synchronous and Poisson gauges. 

We can start by writing the second-order momentum-conservation equation [61] = 0, which gives 


{V 2 i + Uj 2 i) + ^ + Uj 2 i) = -29id> - - di {vj + Lp^) . (A.26) 

For pure growing-mode initial conditions (p' (x (p and vu oc dip, which makes the RHS of this equation the gradient of 
a scalar quantity; thus, the vector part only contains a decaying solution {V2i + u} 2 i) OC a and we can safely assume 
V 2 i = —uJ 2 i (see also Ref. [62] for a discussion on how to obtain the same results using the vorticity vector). 

To proceed one needs the second-order perturbations of the Einstein tensor, 52G^p. these can be found for any 
gauge in Appendix A of Refs. [8, 23]. 

The second-order ‘momentum constraint’ 52G*o = 87rG52T o gives [28] 


u 


' ^ J 


{p + (5i) v\ + V. 


(A.27) 


The pure vector part of this equation can be isolated by first taking its divergence to solve for the combination 
-(- V]/' and then replacing it in the original equation. One obtains [28] 


— = 167rGa^p {vijdiSi — vudjdi) . 

We can further simplify this equation and write: 

V^W2i = + 29^00) , 


with 


p,. ^ 2g^{z)E(z)f{n^) 

^OmHo (1 -|- z) 

where [29, 30] /(flm) = dlnD^/d\na « Hq is the Hubble constant and 

V"0o = ((VVo)" - dApod^d>^po) . 


(A.28) 


(A.29) 


(A.30) 


(A.31) 


For A = 0, the above expression for uj 2 i reduces to Eq. (6.8) of Ref. [42], noting that F{z) = rj in that case. 


c. Tensor perturbations 

For the second-order tensor modes we use the traceless part of the {i — j) Einstein equations, namely Eq. (A. 10). 
To deal with this equation we replace it in the expression for the combination T — Eq. (2.11), together with the 
expression for the vector mode uj 2 i- After some calculations one arrives at [28] 

{x^2j' + 2Hx^/ - = AV^d^d(,R\5] + - d^djR\) + Ad^d^d^diR^^ , (A.32) 

where the traceless tensor 


R\ = d^pdkp - i [Vpf 5i + AnGa^p (^v{vik - 


= 1 + 


3^^0m (1 + 


d podkPo - g (V(/?o) Sf, 


(A.33) 
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Equation (A.32) can be solved by Green’s method, as the corresponding homogeneous equation is the one for linear 
tensor modes, whose analytical solutions are known in the limiting cases ^ 1 or > 1- In the A —> 0 limit one 
recovers the result of Ref. [42], namely 


X2*j(x, i) 


1 




ji (fc^) A 

krj ) ’ 


where ji are spherical Bessel functions of order ^ and 5ij(k) = f (Pye with 


Sij = V^OoSij + didjOo + 2 {didjifoV'^ipo - didkipod^djipo) . 


(A.34) 


(A.35) 


APPENDIX B: DERIVATION OF THE MULTIPOLES 


In this section we will show how to obtain Eq. (3.22). We plug Eq. (3.12) into Eq. (3.15). In Eq. (3.12) the two 
integrals depend on the time derivative of the second-order gravitational potentials (4> -|- 4')', while the contribution 
(first — order)^jg^^g comes from additional second-order corrections. For simplicity let us compute these two parts 
separately in Eq (3.12) 

1 A7^ 1 A7^ 

2~r 

where 


2 T 


[$' + = 


rvo 


dr] 


-y (Oni - r)gmg'iv) - ‘i.gmg'iv) + B[{ri) + ^g{g)g'{g )) 


+ gl,B{g) (V-" (9Vo9*(po) - 


(B.2) 


<=-n(vo-n) 


here B{g) is given in Eq. (3.14) and 


(first - order)^,,tg = 


+ 


rvo 


■ rv 


/ dg 

g{g)g'{g)(pUg,^) + 

/ dg(p 

g"{g)‘Poig,x) 

^ 'Hm 


-dvo 



2^liVrn) 1 


(B.3) 


where we have used Eq. (2.3) into Eq. (3.13). It is important to notice that in Eq. (B.3) (p* = V3(?7*,x = —n(? 7 o — ??*)) 
and the gravitational potential under a given integral must be evaluated along the background geodesics, so that, for 
example. 


rv rn 

/ dfj(p= dfj If 

'Vo dr]o 


r 

= / dfig'{fj)ipoifj,x) 

d r?n 


ic=-h(r;o-j)) Jna ' lx=-n(r;o-7)) 

Correspondingly in Eq. (3.15) we will have 

+ ^'] + [(first - order)^] . 

Let us consider the term in Eq. (B.2). We express it in terms of its Fourier transform 
1 AT2 


(B.4) 


(B.5) 


2 T 


f d^k /ID \ 

.[$' + ^>'] = J J dg ^—^(Oni - l)gmg'ig) - 2gmg'(?7) + B[{g) + 4:g{g)g'{g)j [(p( 


](k) 


- gmBig)K2i'k) 


^—i'k-n{r]o—v) 


(B.6) 
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where iir 2 (k) is defined in Eq. (3.23) and [i^g](k) is the convolution giving the Fourier transform of (/?g(x). Inserting 
Eq. (B.6) into Eq. (3.15) and using the Rayleigh expansion 


Mkx)Y;^ik)YU^) 


(B.7) 


Im 


with the orthonormality of the spherical harmonics we find 

d^k rvo 


/ ^i u rvo 

J dT]je{k{r]o-v)) 


+ £!<!!)+ 

3 9m 9m gt 9m 


- i?(??)K2(k) 




[v^o](k) 

(B.8) 


where we have used Eq. (2.3). 

Let us now consider the term in Eq. (B.3). We Fourier expand the gravitational potential under each integral and 
we use the Rayleigh expansion (B.7). For the first term in Eq. (B.3) the computation is the same as described above, 
while the other terms has a different dependence on n. For instance, in Eq. (B.3) the integral over g"{v) becomes 


- 4 


fVo 

- p 


/ d9 

/ dgip 

9"{9)To{9,^) 

’ Vm 

J Vo 



x=-n(Tjo-v) 


d?ki d^k 2 g"{ri) P ,^ 9 ( 9 ) 


^ ' {2nY {2nY 


'Vo 


d^ki d^k2 r\ g'Yv) 


-4{47:Y Y II / dTj^-^jL^kiim - v)) 

T . A/T. r „ A/f„ J \ ) \ ) J Vm 


^ (27r)3 (2TrY 

LiMi L 2 M 2 ^ \ J \ / -JTIm 

r giv) 

/ drj - jL2{k2{m - 7?))'y?m(kl)(/?m(k2)yLjMi(kl)YL2M2(k2) • 

't] 0 dm 


(B.9) 


and similar expressions hold for all of the three remaining terms in Eq. (B.3). Eq. (B.9) shows that in general AT/T 
does not depend on h only through a single angle k • n, as it happens in linear theory. 

Thus Eq. (B.3) becomes 

/ d^ k 

d9g{9)g'{9)[9^l]{k)e-^'^--^^--^^ 

+ (4^)2 ^ H 

(B.IO) 


d^ki d^k2 
(27r)3 (27r)3 


(pm(kl)(/9m(k2)ALiL2(A:i, A:2)TLjMi(kl)YL2M2(k2) 


where AL^L^{ki, k 2 ) is the function defined in Eq. (3.29). From Eq. (B.IO) it is straightforward to compute the 
corresponding multipoles using Eq. (3.15). We find 


[(first - order)2] = _47r(-i)^ / f driji{k{rio - 9 )) [y)g](k) (k) 

J J-q^ 9m 


9 i 9 ) 9 '{ 9 ) r. 21 


% 

+ Y H j ('27^)3 (-27^)3 


L\M.\ L 2 M 2 
X ^LiL2{kl, k2)YL^Mi(^l)YL2M2{^2) ■ 


(B.ll) 


The sum of Eq. (B.8) and Eq. (B.ll) gives the non-Gaussian part of the multipoles afY, for the late ISW effect written 
in Eq. (3.22). 
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APPENDIX C: CALCULATION OF THE EARLY INTEGRATED SACHS-WOLFE EFFECT AT SECOND 

ORDER 


We consider a flat Universe filled by radiation (p^) and pressureless matter {pm)- In this Appendix we derive the 
equations (3.42) and (3.45). 

The relation (3.42) is obtained through the same steps of Section 3 in Appendix A. We take the traceless part of the 
(i-j)-component of Einstein equations at second-order 520"^j — = 87rG((52T*^ — ^52T\5j). The second-order 

perturbations of the Einstein tensor G^^ can be found for any gauge in Appendix A of Refs. [8, 23]. For the total 
energy momentum tensor -1- one finds 

(I 2 F J = (^P7<527 + ^Pm^2m) + p^{l + W^)v\^v]^ + p^{l + Wm)vlrnV]m , (G.l) 

where = 1/3 and Wm = 0 are the equations of state of radiation and matter respectively, v\^, v\^ their linear 
(scalar) velocities and ^27 = ^ 2 P-ilP-i- Thus we obtain 


(4/ - $) + ^ (Vp)' + 


(4- - $) -b 29V9 j‘P + 4(p5*5jp 


1 

~4 

zn 

p 



2 


p^(l + w^)v\^v]^ 


-H (9*W2j + -b - (^X2i + 27iX2j “ ^^X2i) — 

~ 4“ ^7)^17 4” Pm(X 4- — —Pm(l -b Wm)'l’im 5 


(C.2) 


where p = p-^ + pm and v\^ = v\^vl^. We now apply the operator did^ to get rid of the vector and tensor modes, and 
to solve for the combination (4' — 4)) 


where Q is defined by 


with 


4- - 4> = -4(p2 + Q , 

V2Q= -P-b3A, 

P = Pb 

^ ^ I 1 


(C.3) 


(C.4) 


(C.5) 


where 


and the quantity N is given by 


(C.6) 


= d,d^P] . (C.7) 

In order to get the explicit expression for Q, Eq. (3.43), we need to know the linear velocities and For 
adiabatic initial conditions on large scales (deep in the radiation dominated era) one has that initially (see, 

e.g., Ref. [43]). Moreover it can be shown that on large-scales for adiabatic perturbations the velocities then remain 
equal (see, e.g. Refs. [53, 64]). Therefore we set = uj and from the momentum constraint (O-i) Einstein 

equation we have the expression of the velocities in term of p 


di{p' + Up) = ^*7 + — 

2 \z p ^ p 




where we have used that w.y = 1/3 abd Wm = 0. 
Using Eq. (C.8) in Eq. (C.6) brings 


PP =- 

^ 2TTGa'^p 


1 + 


i £ 2 . 

3 p 


d^if'djp' + H{d^(pdjp' + d^if'djif) + 




Pj 


d''(pdjip 


(C.8) 


(C.9) 
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We now make use of the evolution for the linear growing mode of the gravitational potential for adiabatic perturba¬ 
tions [55] 


ip = -F{r])Ci , 


(C.IO) 


where Ci is the curvature perturbation which is a constant on large scales for adiabatic perturbations, Ci = const. 
Equivalently we can use Eq. (3.32), with The function F{ri) is given in Eq. (5.19). 

Using Eq. (C.9) we thus hnd 


^ [1 + 




1 


P = ^ 1 + ^ ^ 


P J F^ 


Py 


Pj 


F' 


2F- 


.F' 1 


n 


.F' 1 


— +2P—+ - 5+^^ F 


n 2 


P'y 


PJ 


F^ 


V^did^{d^p^djpt), 


d''psfdip<. 


(C.ll) 

(C.12) 


and from Eq. (C.4) 


e = -=li + -l 


-1 


fF' 


\n 


F' 1 


— + 2F— + - 5 + ^ F 


n 2 


P'j 


P 


. (C.13) 


Notice that the second squared bracket is 3/C/lO, where JC is given by Eq. (1.3). Finally by introducing the ratio 
y = Pm!P~i and using d/dlna = H~^d/dr] we can write 


Q = - 


10F2 


QIC, 


(C.14) 


where Q is given by Eq. (3.44). 

We now describe in some details how to obtain the equation (3.45) for the gravitational potential 4'. We start from 
the expression for the curvature perturbation C 2 in the Poisson gauge, see Eq. (2.13), 

C 2 =- 4 /- H^ + AC 2 , (C.15) 

P 

and the (0-0) Einstein equation on large scales 


37i:(4'' -f Tfd-) = ?>'kG52T\ 

In Eq. (C.15) the second-order corrections AC 2 are given by 

AC 2 = 2H^^ + 2^{p' + 2Hip) - 
P P P 


+ i2n^p^ + ipF . 



(C.16) 


(C.17) 


while the second-order perturbations of the Einstein tensor can be found for any gauge in Appendix A of 
Refs. [8, 23]. The (0-0) component of the energy momentum tensor reads 

= -]^52P - Pmvlra - ^P-fVl^ , (C.18) 

where p is the total energy density. We thus express 52P in terms of the gravitational potential through the (0-0) 
Einstein equation, and we insert it in Eq. (C.18) to find 

C 2 = -^ + 24(4-' -f n^) - ^ ^ (C.19) 

p' p' rip' 

where we have dropped the velocities appearing in Eq. (C.17) since they are negligible on large scales. Finally we use 
the relation (C.3) found previously between the two gravitational potentials to eliminate 4> in Eq. (C.19) 


C 2 = -4- + 24(^' + HT) - 2^nQ - + AC 2 . 

p' p' rip' 


(C.20) 
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Switching time variable from d/drj to ’ = d/dlna = H ^d/drj we notice that 


and thus Eq. (C.20) becomes 


- i-5' + (4' + ^') 

2p 


a 


a 


® T 

' = )4+ 

<p" + Q-i^AC2 

VVp \ 

2p 

L 2p J 


(C.21) 


(C.22) 


In the rest of this Appendix we compute the expression (3.46) for the quantity AC 2 defined in Eq. (C.17). This is 
obtained by evaluating AC 2 in the case of radiation plus pressureless matter. 

First of all from the the energy continuity equations p'^ = —AHp^ and = —3Hpm one obtains for the total 
energy density p = p^ + Pm 

^ ^ + 16p^ 

p' H 3pm + ^P-y ' ^ ’ 


and hence 


n— -W - 2H^ = -H^ 6 - 


3pn 


3pm + 


(C.24) 


which is the combination appearing in Eq. (C.17). 

We can now express Sip and Sip' in terms of the graviational potential ip using the total energy continuity equation 
and the (0-0) Einstein equation in the large scale limit. The total energy continuity equation on large scales reads 

Sip' = —3Ti{Sip + Sip) — 3(p -I- p)^ = 0, (C.25) 

where p = + Pm = P-y/3, and dip = Sip-y -I- Sipm, Sip = Sip-y -|- SiPm = Sip-y/3 are the total energy density and 

pressure perturbation respectively. We now take advantage of the adiabaticity of the perturbations 


SlP^ 4 SiPn 

P 7 3 Pm 


to find from the continuity equation 


= -3H 1 + - 


S_i^ 

P' 

Plugging Eqs. (C.24) and (C.27) into Eq. (C.17) yields 
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1 Sp^ \ Sip p' 


3 Ap^ + 3pm J p' U 


AC2 = -n^ 


Sip 
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1-k 


4P7 


Ap^ 3pm 

Finally we use the first-order (0-0) Einstein equation on large scales 


) + - 


p +Hp = --n —, 
2 p 


to get 


AC 2 = 
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(C.27) 

(C.28) 

(C.29) 

(C.30) 


where a dot stands by d/dlna. 

We now use Eq. (3.32), which gives linear growing mode for the gravitational potential, to find 


= J- 
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(C.31) 
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where y = Pm/P-y- With the definition 


we get Eq. (3.46). 


R{y) 


1 + y 
4 + 3y 



4 + 3y; ’ 


(C.32) 
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